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HADRON-PRODUCTION IN HIGH-ENERGY COLLISIONS 



It is a well-known fact that hadrons (for example pions, kaons, protons, neutrons, ■ • ■ , 
and their antiparticles) can be produced in lepton-, hadron- and nucleus- induced reactions 
at sufficiently high incident energies. In such processes, we are witnessing an impressive 
demonstration of energy-conversion into matter, where energy-momentum conservation re- 
quires that the total (c.m.s.) energy of the colliding system must be high enough to creat 
the masses of the produced hadrons and provide them with sufficient kinetic energies. A 
vast amount of data exists for such hadron-formation processes where the multiplicity-, 
transverse-momentum- and rapidity-distributions of the produced hadrons in high-energy 
hadron-hadron, hadron-nucleus and nucleus-nucleus collisions have been accumulated, and 
part of them, especially those at extremely high energies are taken from cosmic-ray experi- 
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ments. 

Much effort has been made to describe these data. The conventional way of doing this is 
to divide such a hadron-production process, in accordance with the currently most popular 
picture for subnuclear structure and subnuclear dynamics, conceptually into three steps: 

In step I, the incident-hadron (a free hadron or a bounded hadron namely that inside 
a nucleus) is pictured as a large swarm of partons (quarks, antiquarks and gluons) where 
everyone of them carries a fraction (known as the Feynman- £ and denoted by xp or simply 
by x, with xp or x between zero and unity) of the longitudinal (that is along the incident 
axis) momentum of this hadron. The distribution of Xp for different kinds of partons (for 
example: the u, d, s valence quarks, the sea-quark-pairs; or the gluons) has been care- 
fully parameterized; and the corresponding parton distributions can be readily found in the 
literature 

In step II, it is envisaged that some of the partons in the projectile-hadron are scattered 
by some of the partons in the target-hadron. They interact with one another according to the 
Feynman rules derived from the given QCD-Lagrangian. Under the assumption that some 
classes of reaction-mechanisms are more important than the others and the perturbative 
methods are applicable, one can evaluate the corresponding Feynman graphs. Details of 
such calculations have been worked out, and can be readily found in the literature 

In step III, the observed final-state hadrons are assumed to be directly related to the 
scattered partons mentioned in Step II. The relations between the scattered partons and the 
observed hadrons are parameterized in terms of "fragmentation functions" . Fragmentation 
functions for all possible kinds of partons have been worked out, and they can also be readily 
found in the literature 

In connection with the above-mentioned three-step-approach, it is of considerable impor- 
tance to recall .J] the following facts: (A) Quarks, antiquarks and gluons have not been, 
and according to QCD and Confinement they can never be, directedly measured. (B) Per- 
turbative methods can be used only when the momentum transfer in the scattering process 
is so large that the corresponding QCD running coupling constant is less than unity; but 
the overwhelming majority of such hadron-production processes are soft in the sense that 
the momentum transfer is relatively low. These facts imply that, it has not been possible, 
and it can never be possible, to have an experimental check for the three steps individually. 

In the present paper, we discuss such hadron-production processes from a different view- 
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point and ask the following questions: Suppose we focus our attention only on the directly 
measurable quantities (such as the rapidity distributions of the produced charged hadrons), 
and try to perform a preconception-free data-analysis where we rely as much as possible on 
the relevant knowledge of Mathematical Statistics [such as the Law of Large Numbers and 
the (Generalized) Central Limit Theorem], shall we be able to extract useful information on 
the reaction mechanisms of such hadron-production processes by examining the experimen- 
tal data? If yes, what can we learn from such information? Can such information be used, 
for example, to make predictions for future experiments, and/or to check the consequences 
of the conventional approach mentioned above? 

Since a considerable part of the available first-hand information on hadron-production 
processes are distributions of few directly measurable quantities (such as the multiplicity 
n c h of charged hadrons; the rapidity y or the pseudorapidity rj and the transverse momen- 
tum p±_ or transverse energy E± of such a hadron), concepts and methods in Probability 
Theory and Mathematical Statistics are expected to play a distinguished role in describ- 
ing/understanding the existing data, and in making predictions for future experiments. 



II. FROM RANDOM WALK TO STABLE DISTRIBUTIONS IN FINANCE AND 

ECONOMICS 

It is well-known that fluctuation studies are of considerable importance in non-equilibrium 
as well as in equilibrium systems; and there are several reasons for this. One of them is that, 
deviations for or fluctuations about the mean values are always present, even when the 
system under consideration is in equilibrium. Another is that fluctuation studies provide a 
natural framework for understanding a large classes of phenomena, among which the best 
known phenomenon is "Brownian motion" or "the random walk" — known through Albert 
Einstein's work starting 1905 Q|. 

It is, however, not very well-known (at least not among the particle physicists) that, five 
years earlier, in 1900, a then young French student named Louis Bachelier has constructed 
a random walk theory for security and commodity markets 3]. The essence of Bachelier's 
theory is the following: Consider the successive differences of two adjacent spot prices 



L B (t,T)=z{t + T)-z(t), 



(1) 
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where z(t + T) and z(t) is the spot price of a stock or that of a piece of commodity at the 
end of a period of time t + T and t respectively; and assume that the set of all Ls(t, T)'s 
are statistically independent, Gaussian (normally) distributed random variables with zero 
mean, and variance proportional to the differencing interval T. This assumption is often 
called jj] the Gaussian hypothesis [sj], and the arguments supporting this hypothesis are 
based on the Central Limit Theorem (cf. Appendix A and the references given therein). 
To be more precise, if the price changes from transaction to transaction are independent, 
identically distributed random variables with finite variance, and if transactions are fairly 
uniformly spaced through time, the Central Limit Theorem will lead us to believe that price 
changes across differencing intervals such as a day, a week, or a month should be Gaussian 
distributed, because they are simple sums of the changes from transaction to transaction. 
Detailed comparisons between this hypothesis and empirical data have been performed by 
a number of authors [4]. Despite the fundamental importance of Bachelier's work Q], it has 
been become obvious that this hypothesis cannot be right as it stands (some details will be 



given in Section IV). In this sense, Mandelbrot's theory |4j, y, |y| can be understood as an 
improvement and a generalization of Bachelier's j^]. 

The main difference between the theory of Bachelier and that of Mandelbrot is the fol- 
lowing. Instead of Gaussian (normally distributed) random variables with zero mean and 
finite (which can be normalized to unity) variance, Mandelbrot made the assertion that the 
variance of the random variable 

L M (t,T) = In z(t + T) -In z(t) (2) 

can behave as if it is infinite; and consistent with this, he suggested that the Gaussian dis- 
tribution should be replaced by a rich class of probability distributions called stable distribu- 
tions among which Gaussian is a limiting case with finite variance. It is known (cf. Appendix 
A and the references given therein) that such distributions allow skewness and heavy tails, 
and have many intriguing mathematical properties. One of the important properties is that, 
instead of the above-mentioned Central Limit Theorem, they satisfy the Generalized Central 
Limit Theorem (cf. Appendix A and the references given therein). Above all, evidences have 
been found that the empirical distributions conform best to the non-Gaussian members of 
this family of distributions. 

It is of considerable importance to keep in mind that the statistical distributions in 
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Mandelbrot's Fractal Approach to Finance and Economics are stable, stationary, and scale- 
invariant (in other words, scaling). Furthermore, such distributions can have one or several 
of the following properties, (a) Repeated instances of sharp discontinuity can combine with 
continuity, (b) Concentration can automatically and unavoidably replace evenness, (c) 
Non-periodic cycles can automatically and unavoidably follow from long-range dependence. 
While some of the mathematical aspects of the above-mentioned concepts are summarized 
in Appendix A, the methods and results of comparing such concepts with experimental data 
will be discussed in Section IV. 



III. RAPIDITY, LOCALITY AND LIGHT-CONE VARIABLES IN 
HIGH-ENERGY HADRON-PRODUCTION PROCESSES 

Quantitative information about the produced hadrons in high-energy collisions exists 
mostly in form of distributions (histograms) of the measured quantities; and one of the few 
distributions that can be, and has been, well measured is the rapidity (y) or pseudorapidity 
(77) distribution of the electrically charged hadrons produced in such processes. 

The rapidity of an observed hadron of mass m, energy E, momentum p = (p\\,p±) where 
p\\ and p± are the momentum in the parallel and that in the transverse direction of the 
collision axis respectively, is defined as 

y= 1 -H^l). (3) 

y 2 y E-p\\' y ' 

This can also be written as y = hi[(E+p\\)/E±\, where E± = M± = (m 2 + |/3j_| 2 ) 1//2 is known 
as the "transverse energy" or "transverse mass". Due to the fact that the overwhelming 
part of the produced hadrons are energetic pions (the mass of which is negligibly small 
compared with the average value of p\\), and that their transverse momenta are exponentially 
distributed with relatively small average \p±\, the corresponding pseudorapidity 

is often used as a good approximation for the rapidity y (y ~ 77), where E in Eq. (3) is 
replaced by p = (p\\ 2 + Ip^f) 1 ^ 2 - One of the advantages of using 77 instead of y is that the 
former can be expressed in terms of the scattering angle (77 = — In tan |), and thus in order 
to determine 77 experimentally, one only needs to measure 6. Furthermore, we can also write 
7/ ~ 77 — ln[(p + p\\)/ \p±\] and treat \p±\ approximately as a constant. 
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Significant fluctuations have been observed in rapidity-distributions in hadron-production 
processes at cosmic-ray j?J , and at accelerator energies 0, |^| . The question we raise here is 
whether such fluctuations can be directly used to gain information on the reaction mech- 
anisms of the above-mentioned processes in general, and information on the space-time 
properties of such reactions in particular. In order to study such properties, it is useful 
to keep in mind that the overwhelming part of the produced hadrons are energetic pions 
which implies that their locations in space-time are mainly concentrated near the light-cone. 
Hence it is not only meaningful, but also useful to discuss the space-time properties of these 
hadrons in terms of light-cone variables (and/or quantities directly related to them) 

x± = -j={t±x\\), (5) 

P± = ^(E± Pll ), (6) 

where (t,x\\,x±) and (E,p\\,p±) are canonical conjugates to each other with the property 
that their 4-product can be written as 

tE - x\\p\\ —x±-p± — x-p + + x + p- - x_l ■ p±. (7) 

Here, we introduce in analogy to rapidity y defined by Eq. (3), a quantity 

which we call locality. This quantity can, again in analogy to y, also be written as I = 
ln[(t — where s± = (t 2 — a;|) 1//2 stands for the "transverse interval" whereby s 2 = 

t 2 — x| — |x_i_| 2 is the square of the space-time "interval" between the world-point x and the 
origin (0,0,0). Furthermore, we can also define a corresponding "pseudolocality" 

A=iln(^% (9) 

2 r + x\\ 

where r in s 2 = t 2 — r 2 stands for r = (x| + |x_l| 2 ) 1//2 . While the dynamics of a system can be 
described either in terms of space-time variables or in terms of their canonical conjugates, 
it is certainly of considerable interest to have a physical picture in space-time. 



Dealing with su 
tation relations 



Dnuclear phenomena, it is clear that we need to take care of the commu- 
between quantities which are canonical conjugates to each other. From 
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[E,t] = i and [pih^ii] = —i, we readily obtain the corresponding commutation relations not 
only between the light-cone variables x± and p±, but also those between I and y as well as 
those between A and 77. While the derivation and the explicit expressions of these relations 
are not very interesting, one of the consequences should be mentioned: By using a very gen- 
eral form of W. Heisenberg's uncertainty principle developed in 1929 by H. P. Robertson |l2j, 
the following can be obtained. If we are interested in the simultaneous specification of two 
observables A and B, the corresponding operators of which obey a non-zero commutation 
relation [A, B] = iC (i is included for convenience; for A = xn and B = p» we have C = H), 
and if we consider a state \ip > which is normalized but otherwise arbitrary (i.e. not neces- 
sarily eigenstate of either A or B), then, the root mean square derivations of A and B as the 
square- roots of the corresponding quantities satisfy /S.A/S.B > ~ |< C >| = 1 1< [A,B] >|, 
which is an exact and precise form of the uncertainty principle. Note that, no matter whether 
C in the commutation relation is an operator or simply a constant, the right-hand side of 
the above equation is always a constant. By applying this to y and I or directly to 77 and A, 
we obtain 

Ar]A\ ~ K, (10) 

where K is a constant (which is independent of the choice of Lorentz-frames) . As we shall 
see in Section IV, the uncertainty principle and its immediate consequences indeed play an 
important role in interpreting the obtained results. 



IV. PROBING STABILITY, STATION ARITY AND SCALING 



As far as the idea and the mathematical basis are concerned, both Bachelier's Gaussian 
hypothesis a and MandehVot, lly po th es is QQQ have to b e considered as m aste r - 8trok es, 
yet what really counts is whether the relevant hypothesis indeed describes the empirical data. 
In this connection, it has been observed that the Gaussian hypothesis [^] can, for example, 
neither explain the typical erratic behavior of the empirical sample second moments, nor 
reproduce the extremely heavy tails of the empirical distributions. On the other hand, Man- 
delbrot's fractal approach to Finance and Economics 0,15,0] is well supported by empirical 
facts. This conclusion is reached after a series of detailed comparisons between empirical 
data and the characteristic features of the Mandelbrot's theory, namely (A) stability (B) 
stationarity and (C) scaling, has been performed j^O, ||. 
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Taken together with the facts mentioned in Section III, these comparisons have led us to 
the following question: Can the concepts and methods introduced to Finance and Economics 
by Bachelier and by Mandelbrot be helpful in studying and understanding fluctuations in 
subnuclear reactions? 

We recall that pseudorapidity i] (or rapidity y) is a continuous variable, and that the 
multiplicity dN/dr)(r)) of charged hadrons at any r\ (within the allowed kinematical range 
Vmin < V < Vmax) can be very well measured. Hence it is not only meaningful but also 
possible to consider for a given event the multiplicity as well as its relative changes in an 
interval A77 of any size [for details, see below especially item (d)] in the allowed kinematical 
range. Based on the facts mentioned above, it seems tempting to introduce, by analogy to 
Mandelbrot's Lm(^,T) shown in Eq. (2), the quantity: 

ax 1 dN , . . , dN . . 
L(?7, A77) = In — (77 + A77) - In — (77), (11) 

and to examine its properties. What we see, as far as the kinematical limits and the exper- 
imental resolution allow, are the following, (a) The L(i], A?y)'s can be any real number, (b) 
The L(r], Ai]) y s are summable, in the sense that the adjacent 77-interval can be combined or 
divided in an arbitrary manner without having any influence on the definition of L(r], At]) 
given in Eq. (11). (c) The set of L(^,Ar])'s can be viewed as approximately statistically 
independent, identically distributed random variables. It follows then from the Generalized 
Central Limit Theorem (cf. Appendix A and the references given therein) that the only 
possible class of limiting distributions should be stable distributions, which reduce to the 
special case Gaussian if the L(rj, Ary)'s have not only a mean, but also a finite variance, (d) 
Since discontinuities are allowed to occur in Mandelbrot's theory, the bin-size can be chosen 
to be so small that even dN/dr]=0 can be included. It should be pointed out, however, that 
dN/drj = does not occur in the JACEE-data as long as the bin-size is not chosen to be 
smaller than Arj = 0.1 which stands for the resolution power of these experiments. (e)The 
variables rj and A77 in L(rj,Arj) are, in contrast to t and T = At in Bachlier's L B (t,T) and 
Mandelbrot's L^(t,T), not part of space-time, but as we have mentioned in Section III, 
the dynamics of the system of produced hadrons should not depend on whether we choose 
space-time or their canonical conjugates as independent variables. 

In the second part of this section, we perform a detailed comparison between our theoret- 
ical expectations and the available data, and show that the empirical distributions obtained 
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from the relevant experimental data are not only stable but also non- Gaussian. Furthermore, 
we present evidence for stationarity and scale-invariance (scaling). 

We consider the cosmic-ray data obtained by JACEE-Collaboration [7] for Si+AgBr 
at 4 TeV/nucleon and those for Ca+C (or O) at 100 TeV/nucleon. These data have 
attracted much attention not only because they show results of hadron-pro duct ion processes 
in nucleus-nucleus collisions obtained at energies much higher than those obtained from 
accelerator-experiments, but also because they exhibit unusually high multiplicities as well 
as significant fluctuations in rapidity-distributions measured in more than 8 to 10 units of 
77. Among the theoretical discussions on these two single-events jjj, the work by Takagi 
1^| and that by Bialas and Peschanski ^| are the most well-known ones. In fact, their 
13, C| two-component conjecture for single-event rapidity-distributions, in particular for 
those published by JACEE-Collabortion |7J has been adopted by most physicists working 
in this field. According to their [3, [3| conjecture, the observed rapidity-distribution of a 
given event can be, and should be, viewed as a superposition of two distinct parts, namely 
a "statistical" part and a "nonstatistical" or "dynamical" part. In Takagi's paper |13j|. 
the "statistical" part is described "by a smooth function with least oscillation" where he 
"assume (s) rather arbitrarily a function of the form ■ • • " . According to Bialas and Peschanski 

n 

|14j , the "statistical" part of such single-event distributions should be determined as follows: 
"If the inclusive rapidity spectrum (averaged over many events) are not available as e.g. 
when one considers a single high-multiplicity event, it is necessary to invoke one's theoretical 
prejudices on the shape of the inclusive spectrum in order to draw definite conclusions." 

In order to study the "nonstatistical" or "dynamical" part, Takagi [l^J suggested that 
"the rapidity density fluctuation of nonstatistical origin may manifest itself as an oscillatory 
pattern • • • " . In his paper, the fluctuations are analyzed by use of power spectrum. To study 
this part, Bialas and Peschanski [3| propose "to study the dependence of factorial moments 
of the rapidity distribution on the size of the resolution" "by adapting the well-known results 
obtained in the investigation of cascading phenomena and turbulent behavior". Based on 
the above-*— two-cogent eooieeture HQ, the random cascade m odel Q 
and the method of factorial moments Q], a large amount of experimental data has been 
analyzed, and many refined models (including scaling models, multifractal models etc.) have 
been developed. An overview of this development together with the obtained results as well 
as detailed references can be found in the review article by De Wolf, Dremin and Kittel . 
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"Understanding" fluctuations in empirical data by conjecturing that the measured distri- 
bution of a given event consists of two or more components (where every component fulfills 
its designed goal) has also been rather popular in Finance and Economics, especially in the 
pre-Mandelbrot era. Based on such a conjecture, many models have been developed and 
very good fits to the data have been produced. But, according to Mandelbrot [^: "A com- 
mon feature of all these approaches, however, is that each new fact necessitates an addition 
to the explanation. Since a new set of parameters is thereby added, I don't doubt that rea- 
sonable curve-fitting is achievable in many cases. " In fact he clearly expresses his doubts 
about the usefulness of such approach in following terms jj, 0] . "This form of symptomatic 
medicine (a new drug for each complaint) could not be the last word!" "In my view, even if 
an accumulation of quick fixes were to yield an adequately fitting patchwork, it would bring 
no understanding." 

We begin our preconception-free data-analysis by noting that the experimental resolution 
in the given histogram is 0.1, and by choosing this to be the smallest 77-interval, Arj. In doing 
so, the 77-range of the Si+AgBr event (hereafter referred to as JACEE1) from rj = —4.0 to 
4.0 is divided into 80 A^-bins, which lead to a set of 79 L(j], Ar/)'s. Similarly, the grange 
of the Ca+C (or that of Ca+O) event (hereafter referred to as JACEE2) from -5.0 to 5.0 is 
divided into 100 A^-bins; hence there are 99 L(j], Ary)'s. In order to probe whether these 
experimental values for L(rj, A7/)'s can indeed be considered as stable random variables which 
satisfy the Generalized Central Limit Theorem, we make use of the definitions and theorems 
summarized in Appendix A. To be more precise, we check this in three steps. 

In Step (i), we calculate the "running sample mean" of the L(r], Ar]) 7 s 



Here, n stands for the ordering number of the A7/-bins where rji = r] min + (i — l)Ar] and the 
ordering begins from 77 = i] min which is -4.0 for JACEE1, and -5.0 for JACEE2 respectively. 
The results are shown in Fig. 1 and Fig. 2 respectively. They seem to have the tendency of 
approaching zero for sufficiently large n in accordance with the Law of Large Numbers (cf. 
Appendix A). 

In Step (ii), we divide each one of the two L(r], Ar/)-sets (JACEE1 and JACEE2) into two 
distinct groups according to the sign of every individual L(r7, A77), and plot separately, in 
accordance with Definition A. 2 given in Appendix A, the right- and the left-tail distributions 




(12) 



i=l 
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P(L > L + ) and P(L < L~) respectively. The result for JACEE1 is shown in Fig. 3, and 
that for JACEE2 is shown in Fig. 4. The purpose of performing the plots are twofold: 
First, we can see whether the distribution exhibits left-right symmetry and thus determine 
the skewness parameter (3 in case the distribution is indeed stable (cf. Theorem A. 3 in 
Appendix A). Second, we can use them as the basis for a detailed stability-test. See Figs. 5 
and 6 for details. 

In Step (iii), we start with the definition of stability 

S m = c m L + 7 m (13) 

and note that in accordance with Definition A.l and Theorem A.l, the necessary and suffi- 
cient condition for the L(r), Ar/)'s to be stable random variables is that for all integers m > 1, 
there exist constants c m > and 7 m such that Eq. (13) is true, where 



s m = YsHvu&v) 



1=1 



= L(n,mAn) 

, dN, . , dN . . 

= in— — (n + mAn) — in— — (n ) 
drj an 

and the L(rji,Arjys are identical copies of L(rj,Arj) and c m = m 1 ^ with < a < 2 and 
7 m is an arbitrary real number. In order to probe the validity of Eq. (13), we choose the 
distributions mentioned in "=" to be tail distributions. As the first step, we calculate and 
plot the right-hand-side of 

c- m \S m - lm )^L (14) 

for the sets of L(rj, A?7)'s obtained from JACEE1 and that of JACEE2 respectively (they 
are of course the same as those shown in Fig. 3 and Fig. 4. The purpose of reproducing 
them here will become clear below). Next, we set m on the left-hand-side of Eq. (14) to 
be 2, 3, 4, 5 and 6 which corresponds to the n- interval mAn respectively (that is m times 
0.1), and determine in every case the most suitable c m and 7 m (m=2, 3, 4, 5 and 6). We 
then use the method described in Step (ii) to calculate the corresponding c^(S m — 7 m ) + 's 
and c^{S m — 7 m )~'s and plot for every m its right-tail distribution P(c~ 1 [5' m — j m ] > 
c m~[Sm — 7m] + ), and its left-tail distribution P(c~ 1 [S' m — j m ] < c^[S m — 7m] - )- It turns out 
that in these two JACEE-events the maximum of m is m = 6, because, beyond this, the 
sample size would not be large enough to have reasonable statistics. Finally, we compare 
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the results obtained through the procedure mentioned above in the same figure (here we 
see why we need the m = 1 case in the same scale). The results for JACEE1 are shown in 
Fig. 5, and those for JACEE2 are shown in Fig. 6. The striking agreement between the 
distributions of the two sides of Eq. (14) shows that the distributions of random variable 
L(j], Arf) J s under consideration are indeed stable distributions. 

It should be pointed out that, with the help of the mathematical tools quoted in Appendix 
A, further information about these stable distributions can be drawn from the results shown 
in Figs. 1-4: In particular, in accordance with the Law of Large Numbers, Figs. 1 and 
2 strongly suggest the existence of zero mean for JACEE1 and JACEE2. Hence it seems 
meaningful to check, whether the corresponding distributions are symmetric with respect to 
the mean value and thus determine the skewness parameter (3. The fact that the right-tail 
distributions and the left-tail distributions are approximately equal as we can explicitly see 
in Figs. 3 and 4, that is P(L > L+) = P(L < L~) for JACEE1 and for JACEE2 respectively, 
shows that (3 = in both cases. 

Next, we probe stationarity, that is, try to find out whether the sets of L(rj, Ar/)'s ob- 
tained from the ^-distributions measured at different times (or time-intervals) have the same 
statistical properties. Since we do not have many sets of //-distribution data available, what 
we can do at present is only to compare the set obtained from JACEE1 with that obtained 
from JACEE2. In Fig. 7, P(L > L + ) and P(L < L~) for these two events are plotted 
together in the same scale. What we see is that, both the right- and left-tail distributions 
obtained from JACEE1 are very much the same as those obtained from JACEE2 — in agree- 
ment with the theoretical expectations. The fact that these two events occurred at different 
times and in reactions at different energies by using different projectiles and targets makes 
the observed similarity particularly striking! 

Due to the role played by the concept of scale-invariance (scaling) in the fractal approach 
to Finance and Economics and its role played in other Sciences including e.g. Condensed 
Matter Physics, it is quite natural to ask: Are there also empirical evidences, or at least 
indications, for scaling in subnuclear reactions? 

One way of checking this is to recall and to apply the method proposed by Mandelbrot 
[f| in analyzing the spot prices of cotton. He performed such tests in two steps: 

In Step 1, he considers the second moment of the daily change of In z(t) with z(t) the 
spot price. Dividing the period under consideration into 30 successive fifty-day samples 
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which can be numbered in chronological order m, he calculates the second moment of every 
sample, and plot them against m (cf. Fig. 1 of Ref. [6]). Here, he sees enormous variability 
in time of the sample second moment, and no tendency of any limiting behavior. 

In Step 2, he evaluates the cumulated absolute-frequency distribution by making use of 
the above-mentioned figure and plot the result on a double logarithmic paper (cf. Fig. 2 of 
Ref. [6]) in order to check whether the obtained points approximately lie on one straight 
line. This is because, mathematically, any quantity N that can be expressed as some power 
of another quantity s, N(s) = s~ T , has the following property. By taking the logarithm on 
both sides, this equation yields InN(s) = —rlns. Hence the scale invariance can be seen 
from the simple fact that the straight line looks the same everywhere, there are no features 
at some scale which make that particular scale stand out; that is, as Ins varies, InN shows 
no kinks or bumps anywhere. What he sees (in Fig. 2 of Ref [6]) is indeed a straight line! 

Following Mandelbrot's method, we divide the 7/-range in JACEE1 and that in JACEE2 
into successive 20-bin equal size samples and give everyone of these samples an ordering 
number m. Due to the fact that, the mean- values of large samples tend to zero (see Fig. 1 
and Fig. 2), the second moment (which corresponds to that shown in Fig. 1 of Mandelbrot's 
paper 6]) is approximately equal to 



which stands for the variance in these two cases. We plot the variance obtained from JACEE1 
in Fig. 8, and that from JACEE2 in Fig. 9. As we can explicitly see, their behavior are 
indeed rather erratic. The corresponding frequency distributions 



are plotted in Figs. 10 and 11 respectively. It seems that, in both cases the data indicate the 
existence of power-law behavior and thus scale-invariance (scaling). It should be pointed 
out, however, since the number of data points in the JACEE-events are much less than those 
for cotton-prices in Mandelbrot's analyses, not only the sample-size in our case have to be 
smaller, but also that the statistics are not as good as those for cotton-prices. In order to 
amend this deficiency, we propose to use the following alternative method. 

Let us consider (instead of the sample variance given in Fig. 1 of Ref. [6] , and those used 




(15) 



i=l 



P(^20 > S 2o) 
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in Figs. 8 and 9 in this paper) the "running sample variance" 

1 n 

S 2 n (A V ) = — 5>fa, A v ) - L n (Ar])} 2 , (16) 
i=i 

where L n (Ar]) stands for the "running sample mean" shown in Eq. (12) and Figs. 1 and 2. 
The range of the ordering number n is n = 2, 3, - • • ,79 for JACEE1, and n = 2, 3, • • • ,99 for 
JACEE2. The "running sample variance" as a function of the ordering number n is shown 
in Fig. 12 and Fig. 13. Also in these figures, we see that the sequential sample moment 
changes rather erratically with respect to n, and does not seem to tend to any limit. That 
is to say, erraticity and non-existence of tendency are nevertheless characteristic properties 
of the plot, although running samples are in general less independent than those obtained 
by dividing the original set into parts (as it is for example the case in Fig.l of Ref. [6] 
and in Figs. 8 and 9 of this paper). In order to see explicitly that running samples do 
not have much influence on the characteristic features of the original set, we also consider 
in the JACEE1 set of 79 (and similarly in the JACEE2 set of 99) Gaussian 

random variables with zero population mean and unit population variance; and we examine 
in particular the "running sample variance" of the set of 79 (and the set of 99) Gaussian 
random variables. What we see (not shown in this paper) is of course distinct differences 

Q 

between the behavior of Gaussian random variables and the that of JACEE-data 7\. A 
more effective way to examine the scaling behavior of such sets of random variables is to 
plot the frequency distributions 

P{S 2 n > si) 

on log — log papers. Such plots for the JACEE-data |Q] are shown in Fig. 14 and Fig. 15 
respectively. The 78 points obtained from JACEE1 and the 98 points obtained from JACEE2 
indeed lie approximately on straight lines in the log — log plots. Having the well-known 
relationship between scaling, power-law behavior and straight-lines on double logrithmatic 
papers in mind, the straight-line structure in the log — log plots and thus the property of 
scale-invariance (scaling) is evident. In order to see whether (and how much if yes) the 
JACEE-data differ from Gaussian random variables, we also use the above-mentioned sets 
of 79 and 99 Gaussian random variables to do the same kind of plots, and show them in 
the corresponding figure, that is in Fig. 14 and Fig. 15 respectively. For such variables, 
the existence of a scale, namely cr 2 = 1 (recall that we are considering Gaussian random 
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variables with zero population mean and unit population variance) can be clearly seen [cf. 
Eqs. (A14) and (A15) in the appendix]. 

A number of conclusions can be drawn from the empirical results plotted in Figs. 14 
and 15, where once again use have been made of the mathematical tools mentioned in 
Appendix A. First of all, the result obtained from such log — log plots has to be considered 
as an effective indicator for scale- invariance (scaling): This is because, in accordance with 
the Law of Large Numbers for sample variance (see Theorem A. 6 in the appendix), there 
is a profound difference between sets of random variables with finite population variance 
and those without finite population variance. While in the former case, the running sample 
variance S 2 should tend to the population variance a 2 [see Eqs. (A14)] for large n which 
plays the role of a scale, there is no such limiting behavior in the latter case. The empirical 
fact that the frequency distributions of S 2 obtained from JACEE1 and JACEE2 exhibit 
power-law behavior shows that the sets of L(r/, Ar^'s obtained from JACEE1 and JACEE2 
do not have any given scale. It can be readily seen from the corresponding log — log plots in 
which the data-points of each JACEE-event lie on a straight line. This behavior is in sharp 
contrast to that of a sample of Gaussian random variables of the same size. Taken together 
with Definition A. 4, A. 5 and Theorem A. 4, A. 6, the above-mentioned observation shows 
that the population variance of the JACEE-sets should be infinite and that the originally 
allowed range for the characteristic exponent a namely < a < 2 should be narrowed down 
to < a < 2. In other words, the set of L(?7,Ar7)'s obtained from JACEE1 and those 
obtained from JACEE2 are indeed non- Gaussian stable random variables. 

Furthermore, the results of stability tests shown in Figs. 5 and 6 can be, and should also 
be, considered as evidence for scale invariance (scaling), because also they show that the 
relevant scale — here the size of 77-intervals — does not play a role. 

In summary, what we have seen in the second part of this section is that, the analysis 
of cosmic-ray data shows that the L(7],A?7)'s defined in Eq. (11) which are the relative 
changes of the multiplicities of charged hadrons between successive rapidity-intervals A77 can 
be considered as mutually independent random variables, satisfying a non-Gaussian stable 
distribution which is stationary and scale-invariant (scaling). The striking similarity between 
the properties of Mandelbrot's LM{t,T) and those of L(t],At]) introduced in this paper 
suggests that the underlying reaction mechanism(s) of the fluctuations in Financial markets 
and those at the subnuclear level of matter are very much the same. To be more specific, we 
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note that everyone of the Ar/'s in L(r), Ar/) is related to its corresponding AA through Eq. 
(10), and the properties of the L(r), K/AA)'s can be readily expressed in terms of the AA's. 
In particular, since AA decreases with increasing Arj, measurements within larger and larger 
values of Arj correspond to measurements within smaller and smaller values of AA. It means, 
in this context, the validity of stability, stationarity and scaling within larger and larger 
Arj intervals implies the validity of such characteristics at smaller and smaller values of AA 
(intervals in pseudolocality near the light-cone in space-time). Having in mind that stability, 
stationarity and scaling are the fundamental characteristics of Mandelbrot's fractal approach 
to Finance, the result of the present empirical analysis should perhaps be considered as 
an indication that concepts and methods used in Nonlinear Dynamics and/or in Complex 
Sciences should be helpful in describing/understanding such hadronization processes. Studies 
along this line are now underway; the results will be reported elsewhere. 



V. CONCLUDING REMARKS 



The result obtained from analyzing the cosmic-ray data for hadron-production, by using 
the concepts and methods introduced by Bachelier and by Mandelbrot to Finance and 
Economics, shows that it is not only possible but also useful to extract information on 
subnuclear reactions directly from empirical distributions of measurable quantities. It shows 
in particular that the fluctuations in rapidity-distributions in the cosmic-ray data have much 
in common with the fluctuations observed in stock market. The striking similarity seems 
to suggest that they are complex phenomena of the same (fractal) nature. The fact that 
non-Gaussian stable distributions which are stationary and scale-invariant (scaling) describe 
the existing data remarkably well calls for further attention. It would be very helpful to 
have comparison with data taken at other energies and/or for other collision processes. 
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APPENDIX A 

This is a mathematical appendix added to the present paper to make it self-contained. 
It contains a brief introduction to stable distributions and a set of definitions and theorems 
(without proofs) taken from monographs and/or textbooks related to this subject jlfil ]. 

The term "stable distributions" stands for a rich class of probability distributions which 
allow skewness and heavy tails, and have a number of other intriguing properties. One of 
them is the lack of closed formulas for densities and distributions for all but a few; such 
exceptional stable distributions are Gaussian, Cauchy and Levy distributions. 

Definition A.l Let X\,X2, ■ ■ ■ ,X n and X be mutually independent random variables 
with a common distribution P which is not concentrated at one point, where n > 1 and S n 
is the sum of these n random variables S n = X± + X 2 + • • ■ + Xn. Then, the distribution 
P is stable in the broad sense if and only if for each n, there exist constants c n > and 7„ 
such that 

S n = c n X + ln . (Al) 

Here, the symbol "=" means equality in distribution, i.e. both expresses obey the same 
probability law. The distribution P is stable in the strict sense if and only if 7 n = for all 
n. 

Theorem A.l The norming constants c n are of the form c n = n 1 ^ with < a < 2, 
where the constant a is called the characteristic exponent of P. 

Theorem A. 2 A random variable X is stable if and only if X — AZ+B, where < a < 2, 
—l<ft<l,A>Q,B&$l and Z is a random variable with characteristic function 

TXOL 

Eexp(iuZ) = exp(— \u\ a [1 — if3t&n — (signu)}) a^l (A2) 

2 

E exp(iuZ) = exp(— \u\ [1 + if3— (signu) In \u\)) a — 1. (A3) 

7T 

Theorem A. 3 A general stable distribution requires four parameters to describe. The 
commonly used set is {a, ft, 7, 5}. Here < a < 2 is known as the index of stability or the 
characteristic exponent; — 1 < ft < 1 the skewness parameter; 7 > the scale parameter; 
and 5 G 3? the location parameter. The three exceptional distributions mentioned at the 
beginning of the appendix, namely 
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(a) Normal or Gaussian distribution, X ~ X(/i, a 2 ) with a density 

f(x) = -^exp(- ^~^ ), -oo<a;<oo, (A4) 
V27TO- 2cH 

(b) Cauchy distribution, X ~ Cauchy^, 5) with a density 

1 7 

= TT1 1\2> ~°° < x < 00 ' ( A5 ) 

7r 7^ + (rr — dy 

(c) Levy distribution, X ~ Levy(j, S) with a density 

/(x) = V^F^ exp( "2T^ ) ' 5< " <00 ' ^ 

are special cases of stable distributions. Also their characteristic functions have precisely 
the form given in Eqs. (A2) or (A3) mentioned in Theorem A. 2. The parameters have 
respectively the following set of values: 

X ~ N(fji, a 2 ) : a = 2, (3 = 0, A = a 2 /2, B = fx; 
X ~ Cauchy('~f, 5) : a — 1, (3 — 0, A — 7, B — 5; 
X ~ Levy( 7 , 5) : a = 1/2, (5 = 1, A = 7, £ = 5. 

Definition A. 2 Let X be a random variable with distribution P(x). Then its right-tail 
distribution T r (x) is 

T r (x) = 1 - P(x) = P(X > x) (A7) 
where x G 3? is any real number. The corresponding left-tail distribution Ti(x) is 

Ti(x) = P(x) = P(X < x). (A8) 

Definition A. 3 A distribution is said to be heavy tailed if its tails are heavier than 
exponential. For a < 2, stable distributions have one tail (when a < 1 and (3 = ±1) or both 
tails (all other cases) that are asymptotically power laws with heavy tails. 

Definition A. 4 Let f(x) be the density function of a stable distribution. The expression 

00 

E(\X\ P ) = J dx\x\ p f(x) (A9) 

—00 

is called the fractional absolute moment of this distribution, where p can be any real number. 
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Theorem A. 4 The fractional absolute moment E(\X\ P ) is finite if and only if < p < a 
where a is the characteristic exponent of this stable distribution (which satisfies < a < 2), 
and that E(\X\) P is infinite for p > a. 

Theorem A. 5 Khinchine's Law of Large Numbers 

Let {Xi(i = 1, 2, • • • , n)} be a sequence of mutually independent random variables with 
a common distribution. If the expectation p = E(Xj) exists, then for every e > 0, the 
probability 

> e} -> for n -> oo. (A10) 
It means in the language of Mathematical Statistics: As n tends to infinity, the expression 

n 

- Xi converges in probability to p\ and this is often written as: 

n i=i 

1 V — "\ p 

— Xi — ► ji for n —> oo. (All) 

i=i 

Definition A. 5 Let {Xi{i = 1, 2, • • • , n)} be a sample of size n taken from the population 
of a random variable X with population mean E(X) = p < oo and population variance 
D(X) = a 2 < oo. Then, the expression 

1 n 

X n = -Vx,, (A12) 

n 

i=i 

which is nothing else but the expression on the left-hand-side of the arrow in (All), and 

1 n 

S n = 7—rJ2^-^) 2 (A13) 

ft -L 

i=i 

are called, respectively, the "running sample mean" and the "running sample variance" of 
X. 

Theorem A. 6 Law of Large Numbers for sample variance 

The sample variance S 2 defined in Eq. (A13) with X n given by (A12) behaves, for 
increasing n, as follows: 

S 2 n a 2 for n -> oo, (A14) 

and for all intergers n, we have: 

E(S 2 n ) = a 2 . (A15) 
Theorem A. 7 Central Limit Theorem 



x 1 + x 2 + --- + x n 



n 
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Let {Xi(i = 1,2, • • • , n)} be a sequence of mutually independent random variables with 
a common distribution. If the expectation /x = E(Xj) and the variance a 2 = Var(Xj) exist, 
then the sum S n — Xi + X 2 H h X„ has the following limiting behavior. For every fixed 

hm P{ ^p^ < X } = ^= I <ftexp(-|). (A16) 

— oo 

This theorem is a special case of the following theorem. 
Theorem A. 8 Generalized Central Limit Theorem 

Let {Xi(i = 1,2, •• • , n)} be a sequence of mutually independent, identically distributed 
random variables. Then, there exist constants a n > 0, b n e 3? and a non-degenerated random 
variable Z with 

a^X, + X 2 + • • • + X n ) - b n -±> Z (A17) 
if and only if Z is stable; in which case a n = rT l l a for < a < 2. 
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FIG. 1: Running sample mean L n (Ai]) of the L(r], A?y)'s [see Eq. (12)] for Ar/ = 0.1 (the experi- 
mental resolution) is plotted as function of the ordering number n. Data are taken from JACEE1 
[3]. Here, r\ m i n = —4.0 is the starting point for the ordering of the 80 A?7-bins. The solid line 
indicates the ideal case in which L n = 0. 
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FIG. 2: Running sample mean L n (Arj) of the L(r], Ar/)'s for Ar/ = 0.1 (the experimental resolution) 
is plotted as function of n. Data are taken from JACEE2 [7]]. Here, rj m i n = —5.0 and there are 100 
Ar/-bins. The solid line indicates the ideal case in which L n = 0. 
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FIG. 3: Tail-distributions of the L(rj, Ar?)'s for JACEE1 pj: The right-tail distribution P(L > L + ) 
is shown as a set of black dots. P(L > \L~\) which is the mirror image of the left-tail distribution 
P(L < L~) is shown as a set of open circles. They are plotted as functions of L + and \L~\ 
respectively. Here, P(L > \L~\) is plotted instead of P(L < L~) in order to avoid the problem 
of taking logarithm of negative values; the values of \L~\ are multiplied by a factor 10 so that 
the two sets of points can be kept away from one another in the figure. This technique used in 
plotting P(L < L~) and comparing it with the corresponding P(L > L + ) in the same figure is 
used throughout the paper whenever it is needed. 
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FIG. 4: Tail-distributions of the L(rj, Ar/)'s for JACEE2 [2j|: Here the right-tail distribution is 
shown as a set of black triangles; while those for the left-tail are shown as open triangles. The 
technique used in plotting P(L < L~) used here are the same as that used in Fig. 3. 
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FIG. 5: Direct stability tests for JACEE1 [2j: Figs. 5(b) to 5(f) correspond respectively to the 
m = 2,3, ••• ,6 cases mentioned in Section IV Step (iii). Here, m stands for the number in 
equal-size-divided samples in the set of L(ry,Ary)'s. This is obviously the same m as that in Eq. 
(14) which expresses the necessary and sufficient condition of stability for the L(rj, Ar/)'s under 
consideration. The black dots and open circles in all these figures are respectively the right-tail 
(r-t) and left-tail (1-t) distributions P{L > L + ) and P{L > \L~\) as functions of L + and \L~\. In 
the figures 5(a) through 5(f) we write simply r-t or 1-t as functions of L + and \L~\ in order to safe 
space. They are to be compared separately with the black stars and the open stars which stand 
respectively for P(c~ 1 [5' m - j m ] > c^[S m - 7 m ] + ) and P(c m 1 [S' m - -y m ] > \c^[S m - 7 m ] _ 
m = 2, 3, • • • ,6. In the figures we write simply r-t or 1-t as functions of = c^~[S m — "f m 



for 
or 



7n 



to safe space. Note that the quality of stability can be judged by the 
degree how good the black stars/open stars overlap with the corresponding black dots/open circles. 
Fig. 5(a) is a reproduction of Fig. 3 in the same scale as Figs. 5(b) to 5(f), so that comparison in 
globe view can be made. 
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(e) (f) 

FIG. 6: Direct stability tests for JACEE2 The notations and the technique in plotting the 
tail-distributions used here are the same as those used in Fig. 5, except that the black and open 
circles should now be replaced by black and open triangles. 
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FIG. 7: Stationarity test: The right-tail (r-t) distribution P{L > L + )and the left-tail (1-t) dis- 
tribution P(L > \L~\) for JACEE1 (abbreviated as Jl in the figure) [J] are plotted in the same 
figure with those for JACEE2 (abbreviated as J 2 in the figure) [jfl. The JACEE1 data are shown 
as black dots and open circles respectively, while the corresponding JACEE2 data are shown as 
black and open triangles. 
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FIG. 8: The variance of samples of the same size (20), Sf (A77 = 0.1), is plotted as function of the 
ordering number m for the L(rj, Arj)'s obtained from JACEE1 jj. Due to the limited number of 
data, m=l, 2, 3 or 4. This figure is the analogue of Fig. 1 of Ref. [6]. 



31 




FIG. 9: Sf (A?7 = 0.1) is plotted against m. Data are taken from JACEE2 p|. Here, m=l, 2, 3, 4 
or 5. Also this figure is the analogue of Fig. 1 of Ref. [6]. 
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FIG. 10: Scaling test for JACEE1 by using the method proposed by Mandelbrot 6]: This is 
the cumulated absolute frequency distribution for the Sfo's shown in Fig. 8. This figure is the 
analogue of Fig. 2 of Ref. [6]. 
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FIG. 11: Scaling test for JACEE2 |jj by using the method proposed by Mandelbrot 6]: This is the 
cumulated absolute frequency distribution for the Sfo's shown in Fig. 9. Also this is the analogue 
of Fig. 2 of Ref. [6]. 




FIG. 12: The running sample variance S^(Ai] = 0.1) for the L(rj, Ary)'s obtained from JACEE1 
is plotted as function of the ordering number n, where n = 2, 3, ■ ■ ■ , 79. 
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FIG. 14: Scaling test for JACEE1 |7j by using the method proposed in the present paper: The 
cumulated absolute frequency distribution for the S^'s (calculated and plotted in Fig. 12) is shown 
as black dots in this figure. For the sake of comparison, a sample of the same size (79 for JACEE1) 
of standard Gaussian random variables (with zero population mean and unity population variance, 
a 2 = 1) is also taken into account, where the running sample variance S 2 's are also calculated, and 
the corresponding cumulated absolute frequency distribution (denoted by open squares) is plotted. 
The solid line stands for the mean of S 2 of the sample of standard Gaussian random variables for 
all n's, that is for E(S 2 ) = a 2 = 1. It also stands for the limiting value of the sample variance, 
which is according to the Law of Large Numbers: S 2 — —* a 2 = 1. 
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FIG. 15: Scaling test for JACEE2 $ by using the method proposed in the present paper: The 
cumulated absolute frequency for the S^'s (calculated and plotted in Fig. 13) is shown as black 
triangles in this figure. The notations used here are the same as those used in Fig. 14. The size of 
the JACEE2 sample and that of the corresponding Gaussian sample is 99. 



